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An effective model that describes the Kondo effect due to a point defect in graphene is developed, taking account of 
the electronic state and the lattice structure of the defect. It is shown that this model can be transformed into a single- 
channel pseudogap Anderson model with a finite chemical potential. On the basis of the numerical renormalization group 
method, it is clarified that the experimentally observed gate-voltage dependence of the Kondo temperature is understood 
in this framework. 
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Recently, the Kondo effect in graphene has attracted much 
attention 1-5 because exotic Kondo effects due to the charac- 
teristic nature of Dirac electrons in graphene 6 are expected. 
It is reported that the Kondo effect is induced by point de- 
fects in graphene. 7 In this experiment, a negative magnetore- 
sistance is observed, which indicates that the Kondo effect is 
the magnetic origin. Since the Kondo effect in usual metals 
is caused by magnetic impurities, this experiment suggests 
that the defects in graphene behave as magnetic impurities. 
The magnetic behavior of the defects in graphene is also sug- 
gested by recent magnetization mesurement of ion-irradiated 
graphene. 8 Although there is no spin at the point defects in the 
usual metals, the magnetic behavior of the defects is known 
in semiconductors, particularly silicon. 9 In graphene, defects 
are introduced by ion irradiation. The most probable defects 
after irradiation are single vacancies, as shown by a scanning 
tunneling microscopy (STM) study. 10 The electronic states of 
point defects in graphene have been studied by density func- 
tional theory (DFT) calculations, 11 ' 12 and the appearance of 
localized spin has been suggested, as in the case of the point 
defect in silicon. 

Although there have been several theoretical studies on the 
Kondo effect in graphene, 1-5 only the adatom cases 1-3 ' 5 and 
the resonant scatterer case 4 have been discussed. No study 
has focused on the localized spin in the point defect itself as 
a cause of the Kondo effect. In addition, experimentally, the 
Kondo temperature depends on the gate voltage. 7 In partic- 
ular, around the charge-neutrality point, the Kondo tempera- 
ture shows a nonzero minimum value and almost symmetric 
dependence for the positive and negative gate voltages. These 
behaviors have not been understood theoretically. 

In this paper, in order to understand the properties of the 
Kondo effect and the related phenomena, we develop an ef- 
fective model in which we take account of the electronic state 
and the lattice structure around the defect. We then analyze 
the model by the numerical renormalization group method. 13 
The relevance of the results to the experimental results is 
discussed. We show that the experimentally observed gate- 
voltage dependence of the Kondo temperature is understood 
in this framework. 

In the case of silicon, the electronic states and the lattice 
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Fig. 1. (Color online) Schematic picture of a point defect (single vacancy) 
in graphene and sp 1 orbitals (circles) around it (<pi, i = 1,2, and 3). -h and 
-h! are transfer integrals between them. 



structures of point defects have been studied in detail. 9 It has 
been known that the electronic state around a point defect can 
be understood by considering molecular orbitals that consist 
of four sp 3 orbitals around the defect. Here, the same treat- 
ment is applied to the point defect in graphene. 

Figure 1 shows the schematic picture of the point defect 
(single vacancy) in graphene. Around the point defect, there 
are three sp 2 orbitals, (p i9 (circles) which do not form cova- 
lent bond, where i indicates one of the three sites (i = 1, 
2, or 3) around the defect. There are also n orbitals on car- 
bon atoms which are not shown in Fig. 1. The hybridization 
between the sp 2 and n orbitals depends on the lattice struc- 
ture of the defect. From STM 10 and transmission electron mi- 
croscope experiments 14 and DFT calculations, 15 two possible 
lattice structures of the point defect are suggested: one is a 
threefold symmetric structure that preserves the original sym- 
metry, and the other is a twofold symmetric structure in which 
one of the C-C distances (between sites 2 and 3 in Fig. 1) is 
shortened. While the former has a planar structure, the latter 
has out-of-plane displacement. When the defect has the pla- 
nar structure, there is no hybridization between the sp 2 and n 
orbitals. However, with out-of-plane displacements, finite hy- 
bridization between the sp 2 orbital, <pi, and the two n orbitals 
at the other sites, j ± i, appears. From the DFT calculations, 15 
it is suggested that the twofold symmetric structure is more 
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stable owing to the Jahn-Teller effect. 

First, we consider an isolated cluster that consists of three 
sp 2 orbitals, as depicted in Fig. 1, as a model for the localized 
state of sp 2 orbitals around the point defect. The tight-binding 
approximation is applied to this cluster. In Fig. 1, -h and -h! 
are transfer integrals between the sp 1 orbitals (h, h f > 0). The 
threefold symmetric structure corresponds to the case with 
h = h', and the twofold symmetric one to h' > h, assum- 
ing that the distance between sites 2 and 3 is shorter than the 
others. 

The energies and the wave functions of this cluster model 
are 
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lff 2 = C+(fi - -^C-((fi 2 + (f3), (1) 



with C+ 
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There are three electrons on 



these sp 2 orbitals. Since the number of electrons is odd, there 
is an unfilled level and nonzero spin appears. Interpreting the 
Kohn-Sham eigenvalues of the DFT calculations 15 in terms of 
these tight-binding energy levels, h'/h is semiquantitatively 
estimated to be h! jh 5. Then, C+ becomes C+ ^ 1, C_ <?c 1, 
and the wave function of the lowest level if/\ becomes the 
bonding state, = {if 2 + $3)1 V2, while the second state 
becomes if/2 = <f\. Thus (f \ becomes the highest occupied 
molecular orbital of this cluster. It should be noted that in a 
divacancy case, four electrons are present in the sp 2 orbitals 
and the localized spin does not appear. 

In the following we consider the structure shown in Fig. 
2, where only one sp 2 orbital (i = 1) is active and it hy- 
bridizes with the 71 orbitals at sites 2 and 3. (The bonding state 
of sp 2 orbitals at sites 2 and 3 is neglected.) We assume the 
hybridization term is given by 



#hyb = Yj [ V ( a L + Ac) d * + h ' C '] ' 



(2) 



<r=U 



where a icr is an annihilation operator of a n electron at site 
i with spin cr =1, I and do- is an annihilation operator of an 
electron on the sp 2 orbital at site i = 1. V is the amplitude 
of hybridization. Taking account of the localized nature of the 
active sp 2 orbital in the defect, we assume that the effective 
Hamiltonian for this orbital becomes 

#def = Y^ 6s p 2 ~ v) d V d <r + Udfadfa, (3) 

cr 

where e sp 2, jj, and U are the energy level, chemical poten- 
tial, and Coulomb interaction on the sp 2 orbital, respectively. 
The Coulomb interaction term is essential for describing the 
Kondo effect. 

The conduction electron states are formed by n orbitals in 
graphene. Ignoring the Coulomb interaction in the n orbitals, 
the conduction electron states are described by the following 
tight-binding model with a nearest-neighbour transfer integral 
K> 0), 6 



(4) 



(ij)cr 



where bjo- is an annihilation operator of a n electron at site j 
of the B sublattice, and N is the total number of n electrons. 
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Fig. 2. (Color online) Schematic picture of the electronic state around the 
defect. Two of the three sp 2 orbitals (i = 2, 3) form a covalent bond. The 
object at site 1 shows the active sp 2 orbital and the circles at sites 2 and 3 
the n orbitals. V is the amplitude of the hybridization between the active sp 2 
orbital and the two n orbitals. 



In eq. (2), only the n orbitals of the A sublattice, a^, appear 
since we have assumed the missing carbon atom to be on the 
B sublattice. Finally, the effective Hamiltonian of graphene 
with the point defect becomes H = H gra + H& Q f + #hyb- 

We apply the numerical renormalization group (NRG) 
method to this effective Hamiltonian to shed light on the elec- 
tronic state at low temperatures. In the usual Kondo problems, 
a magnetic impurity interacts with conduction electrons at the 
same site. However, in the present model, the sp 2 orbital at 
the defect interacts with the n orbitals of sites 2 and 3. Such 
a geometry is treated in the two-impurity Kondo problem. 16 
We follow the method developed for the two-impurity Kondo 
problem. 16 ' 17 

By diagonalizing in k space and retaining low-energy states 
of the conduction electrons, the Hamiltonian of the conduc- 
tion electrons (4) is expressed as 3 ' 6 



gra 



v r 



dA: [s T (k) - p] c\ c kl 



(5) 



with the linear dispersion s T (k) - (-l) T v^|A:|, where vp is 
Fermi velocity, and s(= ±) and r(= 1,2) are valley and band 
indices, respectively. Here, the Dirac points are located at 
K+ = (±4n/3a, 0) in the Brillouin zone corresponding to the 
two valleys, with a being the lattice constant, and k is mea- 
sured from K+ for s = ±. The region of the k integral is lim- 
ited to \k\ < k c , where k c is a cut-off wave number. The anni- 
hilation operator a icr for the A sublattice is expanded in cuo-st 
as 



_ / Qcell y f d r i(k+K s ) ri 

ai(T ~ \ liln) 2 Zi J d * e 



(6) 



where Q ce u is the area of the unit cell of graphene. 

In order to express the hybridization term (2) in a one- 
dimensional representation, we introduce the following op- 



erator: 



.17 



Cka 



(7) 



where S(x) is the Dirac delta function, and i = 2 or 3, r2 = 
R/2, and r 3 = -R/2, with R = (a,0). According to ref. 17 
the following orthonormal states are obtained from a £ 2o- and 

1 



(a £2c r ± a £3(r ), (even, odd) 



(8) 
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where N e (s) and N (s) are given by 

^cell 



N e , (s) 



7lVt 



lT \ J »\T F a 



(9) 



Here, Jo(x) is the zeroth Bessel function. Using these opera- 
tors, the hybridization term (2) can be written as 



#hyb = V Tjf_ D d ^ /2 (*)4rA + h.C 



(10) 



where D = vpk c is a cut-off energy. Note that only one chan- 
nel of the conduction electron (c scre ) hybridizes with the sp 2 
orbital, despite there being four species of conduction elec- 
trons in graphene (s and r). This is because the sp 2 orbital 
hybridizes equally with the n orbitals at sites 2 and 3, as 
shown in Fig. 2. (When the hybridization is not equal, both 
even and odd states are coupled to the sp 2 orbital.) Since 
the dependence of Jo(\s\a/vF) on small \s\alvp is weaker than 
\s\, we approximate Jo(\s\a/vF) as /o(0)(= 1). Assuming that 

a/2 



(9 \ ' 
Q C eii/27rv^J - 1 /D, we obtain 



Hi 



hyb 



^zX/ £|e|i/2 ( c ^ +hx -)- (n) 



Using eqs. (5), (7), and (8), it can be shown that c ecre satisfies 
the commutation relation [//gra, c\ae[ — ( £ ~ v) c lo-e- Thus, we 
assume that the conduction electron term (5) can be written as 



Ho 



ds (s-ii) c\ ae c E[ 



(12) 



Here, we have neglected the other channels of conduction 
electrons. 16, 17 

This model of eqs. (3), (11), and (12) is known as the 
single-channel pseudogap Anderson model 18 with a finite 
chemical potential, where the density of states (DOS) of the 
conduction electron is proportional to \s\. The Kondo prob- 
lem in such pseudogap systems has been studied in detail 18 
for the case with \i - 0. The main feature in this case is 
as follows. In the particle-hole symmetric (2e sp 2 + U = 0) 
case, no Kondo screening occurs, i.e., the localized moment 
is not screened by conduction electrons since the DOS van- 
ishes at id = 0. On the other hand, in the particle-hole asym- 
metric (2e sp 2 + U ^ 0) case, the localized moment is screened 
when the coupling constant V is larger than a critical value 
V c . In real materials, the system is expected to be particle-hole 
asymmetric. In ref. 18, the particle-hole symmetric case and 
the U = oo case were mainly treated, and the case of realistic 
particle-hole asymmetry has not been studied in detail. Hence, 
we focus on the particle-hole asymmetric case with realistic 
parameters. We also study the chemical potential dependence, 
which previously has not been studied in detail. We apply the 
NRG method for the model of eqs. (3), (11), and (12), follow- 
ing the formalism of ref. 18. The NRG method can include the 
effects of the pseudogap DOS and a finite chemical potential. 

In order to clarify the low-temperature states of this model 
in the parameter space of (V, //), we calculate the tempera- 
ture dependence of the entropy of the electron on the sp 2 or- 
bital, Sdef, by the NRG method. 19 The bandwidth of graphene 



is about D ^ 8 eV. We assume that 6, 



-1 eV and 



U ~ 3 eV. 12 Thus we use the parameters e sp 2/D - -0.125, 
and U/D = 0.375, in the following. First, we investigate the 
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Fig. 3. (Color online) Entropy of the electron on the sp 1 orbital, Sd e f, at 
chemical potential \x = as a function of temperature T/D for the hybridiza- 
tion V/D = 0.4, 0.6311, and 0.8. Here, Boltzmann constant k B = 1. 
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Fig. 4. (Color online) Contour plots of Sd e f / In 2 as a function of chemical 
potential fi/D and temperature T/D. 



li ID = case. Figure 3 shows S^f as a function of tempera- 
ture for three values of hybridization V (V/D = 0.4, 0.6311, 
and 0.8). 

We find that the critical value of hybridization is V c = 
0.631 ID. At V = V c , Sdef becomes In 3, which is known as 
the valence-fluctuation fixed point of the pseudogap system. 18 
The Boltzmann constant is set to be ks = 1. For V < V c , 
the entropy remains In 2 at the lowest temperature (~ 10~ 6 D), 
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Fig. 5. Local-moment (LM) region and asymmetric strong-coupling 
(ASC) region in the low-temperature hybridization-chemical potential di- 
agram. The valence-fluctuation (VF) fixed point is located at V c /D = 
0.6311,////) = 0. 
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Fig. 6. (Color online) Kondo temperature as a function of chemical poten- 
tial ii ID for three values of hybridization V/D. Here, e sp i /D = -0.125 and 
U/D = 0.375. 



which indicates an unscreened localized moment. For V > V c , 
Sdef goes to zero, indicating that the localized moment is 
screened. Despite there being no conduction electron DOS at 
[A = 0, the Kondo screening occurs owing to the particle-hole 
asymmetry. 

Figure 4 shows S^f as a function of chemical potential 
ji/D and temperature T/D for three values of hybridization 

V (V/D = 0.4, 0.631 1, and 0.8). In Figs. 4(a) and 4(b), we can 
see that 5def goes to zero at low temperature when ji ± 0. For 

V > V c [Fig. 4(c)], Sdef g° es to zero at low temperatures for 
all values of fi. 

Combining these data, we plot the low-temperature V-\i 
phase diagram in Fig. 5. The whole region is divided into 
the local-moment (LM) region and the asymmetric strong- 
coupling (ASC) region. The valence- fluctuation (VF) fixed 
point is located at V = V c , \i - 0. 

In the experiment, 7 it is observed that the Kondo temper- 
ature depends on the gate voltage. In particular, around the 
charge-neutrality point (fi = 0), the Kondo temperature has 
a finite value, and it shows an almost symmetric dependence 
for the positive and negative gate voltages. 

In order to compare the NRG result with those of the exper- 
iment, we plot the Kondo temperature as a function of \i/D for 
four values of hybridization V/D in Fig. 6. Here, the Kondo 
temperature, T K , is defined by the relation Sfef(T K )/\n2 = 
1 / 2. In the experimental situation, by applying the gate volt- 
age, the chemical potential ji can be varied in the range of 



< 0.2 eV(^ 0.025D). As shown in Fig. 6, T K « 10~ 6 D for 
V < V c at ji = 0. However, T K is finite for V > V c even at 
[A = 0. This behavior reproduces the experimental behavior. 
We find the almost symmetric dependence on the chemical 
potential. Note that this feature is not obtained in the U = oo 
case, 2 where T K becomes totally asymmetric. 

In this study, we have considered only the most stable Jahn- 
Teller distorted structure, i. e., h ± W . However, the threefold 
symmetric structure may be realized under certain conditions. 
In such a case, there is a degeneracy of the molecular orbital s 
at the defect site. It can cause a more exotic Kondo effect, such 
as the three-impurity Kondo effect. 20 This remains as a future 
problem. Also, the dynamical effects of the lattice degrees of 
freedom at the defect site have been neglected. Actually, the 
dynamical switching of the covalent bond (dynamical Jahn- 
Teller effect) in the point defect has been discussed. 15 ' 21 We 
speculate that this effect enhances the electron scattering at 
the defect site and hence enhances the Kondo effect. 22 

In conclusion, a point defect in graphene has been investi- 
gated as a cause of the Kondo effect. Taking account of the 
electronic state and the lattice structure of the defect, an ef- 
fective model has been developed. This model is transformed 
into the single-channel pseudogap Anderson model with a fi- 
nite chemical potential. By numerical renormalization group 
calculations, the low-temperature phase diagram of this model 
is determined in the whole parameter regions of hybridization 
and chemical potential. It is shown that the experimentally 
observed gate- voltage dependence of the Kondo temperature 
can be reproduced in this framework. 
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